Gamma Ray Bursts as standard candles to 
constrain the cosmological parameters 



G. Ghirlanda, G. Ghisellini, C. Firmani 

Osservatorio Astronomico Brera, Via Bianchi 46, 1-23807 Merate (LC), IT 
E-mail: giancarlo . ghirlandaObrera. inaf . it 

Abstract. Gamma Ray Bursts (GRBs) are among the most powerful sources in the 
Universe: they emit up to 10 erg in the hard X-ray band in few tens of seconds. The 
cosmological origin of GRBs has been confirmed by several spectroscopic measurements 
of their redshifts, distributed in the range z £ (0.1,6.3). These two properties make 
GRBs very appealing to investigate the far Universe. Indeed, they can be used to 
constrain the geometry of the present day universe and the nature and evolution of 
Dark Energy by testing the cosmological models in a redshift range hardly achievable 
by other cosmological probes. Moreover, the use of GRBs as cosmological tools could 
unveil the ionization history of the universe, the IGM properties and the formation of 
massive stars in the early universe. The energetics implied by the observed fluences and 
redshifts span at least four orders of magnitudes. Therefore, at first sight, GRBs are 
all but standard candles. But there are correlations among some observed quantities 
which allows us to know the total energy or the peak luminosity emitted by a specific 
burst with a great accuracy. Through these correlations, GRBs becomes "known" 
candles, and then they become tools to constrain the cosmological parameters. One of 
these correlation is between the rest frame peak spectral energy E pca k and the total 
energy emitted in 7-rays Ej, properly corrected for the collimation factor. Another 
correlation, discovered very recently, relates the total GRB luminosity Lj so , its peak 
spectral energy -Speak and a characteristic timescale J0.45, related to the variability of 
the prompt emission. It is based only on prompt emission properties, it is completely 
phcnomenological, model independent and assumption-free. These correlations have 
been already used to find constraints on and Vt\, which are found to be consistent 
with the concordance model. The present limited sample of bursts and the lack of 
low redshift events, necessary to calibrate the correlations used to standardize GRBs 
energetics, makes the cosmological constraints obtained with GRBs still large compared 
to those obtained with other cosmological probes (e.g. SNIa or CMB). However, the 
newly born field of GRB-cosmology is very promising for the future. 
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1. Introduction 

The present-day Universe and its past evolution are described in terms of a series of 
free parameters which combined together make up the cosmological model. One of the 
key objectives of observational cosmology is to exploit the power of combining different 
observations in order to test different cosmological models and find out the "best" 
one. The standard model of cosmology can have up to about 20 parameters needed 
to describe the background space-time, the matter content and the spectrum of the 
metric perturbations (e.g. PUI2] )• However, a "minimal" subset of 7 parameters can be 
used to construct a successful cosmological model. In addition to the set of parameters 
describing the global geometry and dynamics of the Universe (in terms of its curvature 
and expansion rate - i.e. Qm, and h), it became of greatest interest also the 

description of how the structure were formed from the basic constituents (described in 
terms of the perturbations amplitude A), the ionization history of the universe since 
the era of decoupling (r) and the description of how the galaxies trace the dark matter 
distribution on the largest scales (described in terms of the bias parameter b). 

The goal of observational cosmology is to use astronomical objects and observations 
to derive the cosmological parameters (i.e. cosmological test - e.g j3l Ej ) once a 
cosmological model has been defined by selecting a set of parameters (i.e. the model 
selection problem - e.g. [JJE])- 

The Hubble constant h = 0.72 ± 0.1 has been recently measured by the HST key 
project through the calibration of primary (Cepheid variables) and secondary (SN type 
I, SN type II, Galaxies) distance indicators within 600 Mpc. This measure has been 
confirmed by the CMB WMAP data and by large scale structure measures (jH|,[I|). 

Supernovae la, through the luminosity distance test (see below), have been used to 
constrain (combined with the evidence of a flat universe from CMB data) Qm — 0.3 and 
£7a — 0.7 (e.g. [HI HH1 HU IE])- One of the greatest breakthrough of the cosmological 
use of SNIa is the evidence of the recent re-acceleration of the universe (e.g. [TB^ ITU ITS] 
see also ^Ej for a review), interpreted as the effect of a still obscure form of energy with 
negative pressure. 

The CMB primary anisotropies bear the imprint of the physical conditions at the 
epoch of matter-radiation decoupling and of several effects such as the time-varying 
gravitational potential of structures along its propagation path (the ISW effect), the 
gravitational lensing and the scattering by the homogeneous ionized gas and by the 
collapsed-ionized gas (the SZ effect). Among the greatest breakthroughs of the CMB 
data analysis ([Zj) it is worth mentioning that (i) the universe is flat, at the 1% level of 
accuracy, (ii) r ~ 0.09 ± 0.03, (iii) the power spectrum of the initial perturbations has 
the form of a scale invariant powerlaw, (iv) a static dark energy model is preferred. In 
particular the combination of the three year WMAP data with the Supernova Legacy 
data ([12]) yields a significant constrain on the dark energy equation of state, i.e. 
w = —0.97 ± 0.02, favouring the ACDM cosmological model. 

The study of the power spectrum of large galaxy surveys allowed to derive the 
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baryon fraction Q^/Qm = 0.185 ± 0.046 (|17j) or, independently from the CMB priors, 
Qm = 0.273 ± 0.025 ([HI)- Galaxy clustering analysis has also been used to put 
constraints, within the CDM model, on the neutrino mass (e.g. [19J). On the other 
hand Galaxy clusters, combined with constraints on the baryon density from primordial 
nucleosynthesis, have been used to constrain 0« and h Q20J). Finally, strong (by 
clusters potentials) and weak gravitational lensing have been used to measure the mass 
power spectrum and derive constraints on a 8 and Qm (e.g. |2*T]). 

Gamma Ray Bursts are very promising for cosmology for, at least, two reasons: (i) 
they cover a very wide redshift range (see Fig. Q solid black line) presently extending 
up to z = 6.29 ([22], see [B]) and (ii) GRBs are detected by space instruments in the 
7-ray band at energies >10 keV, i.e. their detection is free from the typical limitations 
due to dust extinction in the optical band. 

The prospects of using GRBs as a new cosmological tool are exciting: 

• similarly to QSO, GRBs might contribute to study the distribution and the 
properties of matter observed in absorption along the line of sight to these distant 
powerful sources (e.g. [23] , [21] , [213 , [2H] ) • Moreover, given the high luminosity of the 
early afterglow and the GRB transient nature, the absence of "proximity effects" 
(typical of QSO) makes GRBs powerful sources to study the metal enrichment and 
ISM properties of their own hosts ((27j.|5Sj. P] see also [H]); 

• if GRBs correspond to the death of very massive stars, they represent a unique 
probe of the initial mass function and of the star formation of massive stars (|31j) 
at very high redshifts; 

• GRBs can be potentially detected at any redshift: the spectroscopic study of GRBs 
will help to define the properties of the IGM beyond the present QSO limits (z ~ 6) 
and, consequently to study the epoch(s) of re-ionization (e.g. |SSj) 

• the correlations between GRB spectral properties and collimation corrected 
energetics ([M],[II2]), among prompt observables (jSH]) and prompt and afterglow 
observables ([39J) have been shown to be powerful tools that "standardize" GRB 
energetics which can therefore be used to constrain the universe dark matter and 
and energy content (f2 M ,^A) and the nature of dark energy ( [3"4] . |4T] . [4*2"] . [313] . |4*3] . 

W)- 

GRBs are not alternative to SN la or other cosmological probes. On the contrary, 
they are complementary to them, because of their different redshift distribution and 
because any evolutionary properties would likely to act differently on GRBs than on SN 
la or other probes and because the joint use of different cosmological probes is the key 
to break the degeneracies in the found values of the cosmological parameters. 

1.1. The standard candle test 

The classical methods used to test the cosmological models are: (i) the luminosity 
distance test (mainly applied to SN la), (ii) the angular-size distance test (whose modern 
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application consists in the study of the CMB anisotropies) and (hi) the volume test 
based on galaxy number count. Besides there are cosmological tests aimed at testing 
structures' formation models such as the study of the power spectrum of luminous 
matter or of the amplitude of the present-day mass fluctuations. 

For an object with known and fixed luminosity L, not evolving with cosmic time, 
and measured flux F, we can define its luminosity distance Dl = (L/AttF) 1 ^ 2 which 
is related to the radial coordinate r of the Friedman-Robertson-Walker metric by 
Dl = r/a{r) = r (1 + z) (where a(r) is the scale factor as a function of the comoving 
time, t = tH ). Therefore, Dl depends upon the expansion history and curvature 
of the universe through the radial coordinate r. By measuring the flux of "standard 
candles" as a function of redshift, F(z), we can perform the classical luminosity distance 
test by comparing Dl, obtained from the flux measure, with Dl(p) predicted by the 
cosmological model, where p is a set of cosmological parameters (e.g. Qm, and h). 

This test has been widely used with SNIa (see [35 J for a recent review). The high 
peak luminosity (i.e. ~ 1O 1O L ) of SNIa makes them detectable up z >1. However, 
SNIa, strictly-speaking, are not standard candles: their absolute peak magnitude varies 
by ~ 0.5 mag which corresponds to a variation of 50%-60% in luminosity. In the 
early 1990s (j3H|, jUj, |H],[39J) it was discovered a tight correlation between the peak 
luminosity of SNIa and the rate at which their luminosity decays in the post-peak 
phase. This is known as the "stretching-luminosity correlation" and in its simplest 
version^ is M B ~ 0.8(Am 15 — 1.1) — 19.5, where M B is the B-band absolute magnitude 
and Amis represents the decrease of the magnitude from the time of the peak and 15 
days later. The application of this correlation reduces the luminosity spread of SNIa to 
within 20% and makes them usable as cosmological tools ( 48J). 

The situation with GRBs is remarkably similar. At first glance the wide dispersion 
of their energetics prevented their use as cosmological probes even when accounting 
for their "jetted" geometry (Sec. 2). Several intrinsic correlations between temporal or 
spectral properties and GRB isotropic energetics and luminosities (Sec. 3) have been 
regarded as a possibility to make GRBs cosmological tools. However, only by accounting 
for a third variable (which in the standard model measures the GRB jet opening angle - 
Sec. 4), GRBs became a new class of "standard candles" (or, better, "known" candles). 
The constraints on the cosmological parameters obtained by the luminosity distance test 
applied to GRBs are less severe than what obtained with SNIa (Sec. 5). This is mainly 
due to the presently still limited (i.e. only 19) number of GRBs with well determined 
prompt and afterglow properties that can be used as standard candles (Sec. 9). 

2. GRB energetics and luminosities 

The cumulative redshift distribution of long GRBs is compared in Fig. ^ with that 
of 156 SNIa We also show the different redshift distributions obtained with 

| see also |5()| for and alternative definition of the stretching-luminosity correlation based on multi- 
color modeling of SNIa lightcurves. 
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Figure 1. Normalized redshift distribution of 79 long duration GRBs updated to Jan 
2006 (solid-black histogram) . Also shown are the distributions separated according to 
three different satellites that detected most of these burst (BeppoSAX [15 GRBs, dot- 
dashed blue line], Hete-II [17 GRBs, dashed-green line] and Swift [29 GRBs, dotted- 
red line]). The redshift distribution of the "Gold" SNIa sample (156 objects) of 
is also shown (solid-blue line). The insert shows the sky distribution, in galactic 
coordinates, of the present sample of GRBs with spectroscopically measured redshifts 



(see http://www.mpe.mpg.de~ jcg/grbgen.html for reference) 



the GRBs detected by the three different satellites Swift, BeppoSAX. and Hete-II. The 
better sensitivity and faster accurate afterglow localization of Swift ([HI]), compared to 
BeppoSAX. and Hete-II, might account for the detection, on average, of fainter X-ray 
and Optical afterglows and, therefore, for systematically larger redshifts (|S2])- It should 
also be noted that due to its soft energy band (15-150 keV), Swift might better detect 
and localize soft/dim long bursts. The insert of Fig. ^shows that the sky distribution of 
the population of bursts with measured z is consistent with that of long GRBs detected 
by BATSE (see e.g. 

2.1. Isotropic energy and luminosity 

The energy (E) and luminosity (L) of GRBs with measured redshifts can be estimated 
through the burst observed fluence (S, i.e. the flux integrated over the burst duration) 
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and peak flux (P). If GRBs emit isotropically, the energy radiated during their prompt 
phase is E iso = AttD 2 l S/ (1 + z) [where the term (1+z) accounts for the cosmological time 
dilation effect] and the isotropic luminosity is L iso = AnD 2 L P. Fig. El (upper panel) shows 
the distribution of E- lso (orange-filled histogram) obtained with the most updated sample 
of 44 GRBs (jSSj) with measured redshifts and measured spectral properties (from which 
the bolometric fluence in the 1-10 4 keV rest frame energy band can be computed). In 
the same plot we also show the distribution of L iso (blue-hatched histogram) which has 
been obtained by updating the compilation of [HE] with the most recent events. Note 
that the bolometric luminosity is often computed by combining the peak flux (relative 
to the peak of the GRB lightcurve) with the spectral data derived from the analysis 
of the time integrated spectrum. As discussed in jSH] this is strictly correct only if the 
GRB spectrum does not evolve in time during the burst, contrary to what is typically 
observed (|H7j). 

If the GRB energy and luminosity distributions of Fig. |21 (upper panel) are modeled 
with gaussian functions (solid-red line and solid-green line, respectively), we find that 
the average isotropic energy of GRBs is (Log(E iso )) = 53.03 ± 0.89 while the average 
luminosity is {Log(L iso }) = 52.23 ±0.10 (la uncertainty). Clearly, this large dispersion 
of E iso and L iso prevents the application of the luminosity distance test to GRBs. 

2.2. Collimation corrected energy and luminosity 

A considerable reduction of the dispersion of GRBs energetics has been found by [SB] 
(later confirmed by when they are corrected for the collimated geometry of these 
sources. 

Theoretical considerations on the extreme GRB energetics under the hypothesis of 
isotropic emission (|60j.|61j) led to think that, similarly to other sources, also GRBs 
might be characterized by a jet. In the standard scenario, the presence of a jet ([62J, 
[EB] i |EI] affects the afterglow lightcurve which should present an achromatic break few 
days after the burst (see Sec. 4 for more details). 

Indeed, the observation of the afterglow lightcurves allows us to estimate the 
jet opening angle 8j Ct from which the collimation factor can be computed, i.e. / = 
(1 — cos ^jet). This geometric correction factor, applied to the isotropic energies (|5H]) led 
to a considerable reduction of the GRB energetics and of their dispersion. Fig. |21 (lower 
panel) shows the distributions of the collimation corrected energy and peak luminosity 
(Ky = -E; so • / and L 7 = L iso • /, respectively) for the updated sample of bursts with 
measured 6- ieXj . When compared with the parent distributions of the isotropic quantities 
(upper panel of Fig- we note that, also accounting for the collimation factor, both 
and L 7 are spread over ~3 orders of magnitudes. 

However, given the slightly lower dispersion of the collimation corrected energies 
and luminosities (bottom panel of Fig. |2J) with respect to their isotropic equivalent 
(upper panel of Fig. EJ), it is worth to verify if GRBs can be effectively used as standard 
candles. The simplest test consists in building the Hubble diagram (see also [55j.|40]). 
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Figure 2. Upper panel: distributions of the isotropic equivalent energy £"i SO 
(orange-filled histogram) and of the isotropic equivalent luminosity L lso (blue-hatched 
histogram). The sample comprises the 44 long GRBs with measured redshifts and 
well determined spectral properties. The values of -Ei so are taken from I55| while 
Li SO is from (updated with the most recent bursts). The solid lines represent 
the Gaussian fit to these distributions. The two low-luminosity bursts (980425 and 
031203) are also shown (see jSJ, jSl])- Lower panel: distributions of the collimation 
corrected energy (orange-filled histogram) and luminosity (blue-hatched histogram). 
The samples comprises only the 19 GRBs with firm measurements of the jet break time 
from which the jet opening angle could be computed (18 bursts from |35| updated with 
GRB 051022, see g3]). 




E(1 -cos0) 
L(1 -cos5) 



By assuming the central value of the collimation corrected energy, i.e. E 1 = 8.8 x 10 
erg, for all the GRBs with measured z and #j et , it is possible to derive their luminosity 
distance as if they were strictly speaking standard candles, i.e. characterized by a unique 
luminosity. By solving numerically the equation E 1 = 8.8 x 10 50 = E [so ■ [1 — cos0j et ] 
(where also 0j et depends on D L through E iso , see Sec. 4) we derive the luminosity distance 
which is plotted versus redshift in Fig. El Note that the scatter of the data points is larger 
than the separation of different cosmological models (lines in Fig. EJ). This simple tests 
shows that, even if the collimation correction reduces the dispersion of GRB energies, 
they cannot still be used as "standard candles" for the luminosity distance test. 
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Figure 3. Hubble diagram for GRBs with measured redshifts and jet break time 
assuming a unique value of the collimation corrected energy E 1 = 8.8 x 10 50 erg, i.e. 
the central value of its distribution (blue-hatched histogram of Fig|31 lower panel). 
The luminosity distance is obtained by solving the equation _E 7 =const=i?i so (l — cos#). 

3. The intrinsic correlations of GRBs 

For GRBs with known redshifts we can study their rest frame properties. Although 
still based on a limited number (few tens) of events, this analysis led to the discovery 
of several correlations involving GRBs rest frame properties. In the following sections 
these correlations are summarized. 

3.1. The Lag-Luminosity correlation (r-L iso ) 

The analysis of the light curves of GRBs observed by BATSE in 4 broad energy ranges 
(i.e. roughly 25-50, 50-100, 100-300 and >300 keV), led to the discovery of spectral lags: 
the emission in the higher energy bands precedes that in the lower energy bands ( [65J ) . 
The time lags typically range between 0.01 and 0.5 sec (even few seconds lags have 
been observed - jHHl) and there is no evidence of any trend, within multipeaked GRBs, 
between the lags of the initial and the latest peaks ([HZ|)- It has been proposed that 
the lags are a consequence of the spectral evolution (e.g.[p f Hj). typically observed in 
GRBs ([69 ), and they have been interpreted as due to radiative cooling effects Q70J). 
Alternative interpretations invoke geometric (i.e. viewing angle) and hydrodynamic 
effects ([H], [121) w hhin the standard GRB model. 

In particular, the analysis of the temporal properties of GRBs with known 
redshifts revealed a tight correlation between their spectral lags (r) and the luminosity 
(L iso ) Q73J): more luminous events are also characterized by shorter spectral lags as 
represented in FigEI (left panel) where the original correlation r oc L~ o 8 , found with 6 
GRBs, is reported (see also Tab.l). 
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Table 1. GRB intrinsic correlations involving isotropic quantities. N: number of 
GRBs used to find the correlation. The rank correlation coefficient r s and its chance 
probability P are given. The best fit parameters of the powerlaw model (normalization 
q and slope a), with their la uncertainty, are also reported together with the x 2 (dof). 



Isotropic Correlations 


N 




P 




a 


x 2 (dof) 


Lag-Luminosity 


6 


-0.82 


0.06 


39.09±7.2 


-0.77±0.13 




Log(r Ug ) = q + aLog(L iso ) 














Variability-luminosity 


31 


0.61 


8xl0~ 4 


-27.3±1.3 


0.5±0.02 


1252(29) 


Log(V) = q + aLog(L iso ) 














Peak energy-Isotropic energy 


43 


0.89 


lO" 15 


-27.64±0.65 


0.56±0.01 


493(41) 


Log(E pca _ k ) = q + aLog(E iso ) 














Peak energy-Isotropic luminosity 


10 


0.86 


io- n 


-27.02±0.56 


0.50±0.01 


454(38) 



Log(E pcak ) = q + aLog(L iso ) 



Possible interpretations of the Lag-Luminosity correlation include the effect of 
spectral softening of GRB spectra ([Hj, [IS]) during the prompt due to radiative cooling 
([IE] j |ZZ|) or a kinematic origin due to the variation of the line-of-sight velocity in 
different GRBs ([78j) or to the viewing angle of the jet ([79*]). 

Moreover, the r — Li SO correlation has been used as a pseudo-redshift indicator to 
estimate z for a large population of GRBs (jHOl) and also to study the GRB population 
properties (i.e. jet opening angle, luminosity function and redshift distribution) within 
a unifying picture (|66j). 

3.2. The Variability-Luminosity correlation (V -L lSQ ) 

The light curves of GRBs show several characteristic timescales. Since their discovery, 
it was recognized that the 7-ray emission can vary by several orders of magnitudes (i.e. 
from the peak to the background level) on millisecond (or even lower (jH2j) timescales 
(e.g [HI])- Also the afterglow emission presents some variability on timescales of a few 
seconds (e.g. [EH], [EI]). The analysis of large samples of bursts also showed the existence 
of a correlation between the GRB observer frame intensity and its variability (|81j). A 
short timescale variability during the prompt emission of GRBs has been considered as a 
strong argument against the external shock model, favouring, instead, an internal shock 
origin for the burst high energy emission (e.g. [EH])- in fact, the rapid variation of the 
emission was interpreted as a signature of a discontinuous and rapidly varying activity of 
the inner engine that drives the burst ([EH])- This is hardly produced by a decelerating 
fireball in the external medium. However, alternative scenarios propose an external 
origin of the observed variability as due to the shock formation by the interaction of the 
relativistically expanding fireball and variable size ISM clouds (|87|). 

Fenimore & Ramirez-Ruiz (jHH]) and Reichart et al. ( 89J ) found a correlation 
between GRB luminosities (L; so ) and their variability (V): more luminous bursts have a 
more variable light curve. The V — L iso correlation has been recently updated ([MI]) with 
a sample of 31 GRBs with measured redshifts. This correlation has also been tested 
(|91j) with a large sample of 551 GRBs with only a pseudo-redshift estimate (from the 
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Figure 4. Left: Lag-Luminosity correlation (data from |55]) for 6 GRBs with 
measured redshift. The best fit (solid line) is also shown. Right: Variability- 
luminosity correlation for the sample of 31 GRB of 90 . The violet-shaded region 
represents the 3cx scatter of the data points around the best correlation (solid line). 
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Figure 5. Left: Peak energy-isotropic energy correlation. Data are from |55| . The 
orange shaded region represents the 3tr scatter of the data points around the correlation 
(solid line). Right: Peak energy-isotropic luminosity correlation. Data (40 GRBs) 
are from updated with the most recent bursts. The blue shaded region represents 
the 3er scatter of the data points around the correlation (solid line) . 



lag-luminosity correlation - [HO])- An even tighter correlation (i.e. with a reduction 
of a factor 3 of its scatter) has been derived (|S2l) by slightly modifying the definition 
of the variability first proposed by jEHj- In FigJUwe show the data points (from [90] ) 
which show the existence of a statistically significant correlation (with rank correlation 
coefficient r s = —0.6 and chance correlation probability P = 8 x 10~ 4 ). We fitted this 
correlation by accounting for the errors on both coordinates and the best fit coefficients 
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are reported in Tab.l. The scatter of the data points around this correlation (computed 
perpendicular to the correlation itself) can be modeled with a gaussian with a = 0.24. 
As also shown in Fig0] (right panel) such a large scatter is responsible for a statistically 
poor fit, as shown by the resulting x 2 =1252/29. However, as discussed in j^H] when 
the sample variance (see also |HH] , [HI] ) is taken into account (cr; 09 y = 0.20 ± 0.04) the 
correlation is L oc \z 3 - 4 (+°- 9 >-°- 6 ) with a reduced \\ ~ 1- 

3.3. The spectral peak energy-is otropic energy correlation (E pea k-Ei SO ) 

It has been shown ([HE!) that GRBs with highly variable light curves have larger vF v 
spectral peak energies. The existence of a V-E pea ^ correlation and of the variability- 
luminosity correlation (V-Li SO ) implies that the rest frame GRB peak energy Speak is 
correlated with the intrinsic luminosity of the burst. Amati et al. ( 96 ) analyzed the 
spectra of 12 BeppoSAX GRBs with spectroscopically measured redshifts and found 
that the isotropic-equivalent energy E iso , emitted during the prompt phase, is correlated 
with the rest-frame peak energy of the 7-ray spectrum E peilk = -E°^ k (l + z). Such a 
correlation was later confirmed with GRBs detected by BATSE, Hete-II and the IPN 
satellites and extended with X-ray Flashes (XRF) towards the low end of the E peak 
distribution ( , [M] , [EE] ) • Recently the E pea k-E- lso correlation has been updated with a 
sample of 43 GRBs (comprising also 2 XRF) with firm estimates of z and of the spectral 
properties ([EE])- In Fig. El we show the E peak -E iso correlation with the most updated 
sample of GRBs reported in [EE]- The best fit correlation parameters are reported in 
Tab.l. The scatter of the data points around the E pea k-E iso correlation can be modeled 
with a Gaussian of a ~ 0.15. 

The theoretical interpretations, proposed so far, of the E peak -E- lso correlation ascribe 
it to geometrical effects due to the jet viewing angle with respect to a ring-shaped 
emission region ( |HH] , ISO] ) or with respect to a multiple sub-jet model structure which 
also accounts for the extension of the above correlation to the X-ray Rich (XRR) and 
XRF classes (400, [101 J. An alternative explanation of the E pea ^-E- lso correlation is 
related to the dissipative mechanism responsible for the prompt emission Q102J): if 
the peak spectral energy is interpreted as the fireball photospheric thermal emission 
comptonized by a dissipation mechanism (e.g. magnetic reconnection or internal shock) 
taking place below the transparency radius, the observed correlation can be reproduced. 

3. 4- The peak spectral energy-is otropic luminosity correlation (E pea k-L- lso ) 

A correlation between -E^k and the isotropic luminosity L iso (i.e. -Speak oc L^) has been 
discovered ( |109j ) with a sample of 16 GRBs. A larger sample of 25 GRBs confirmed 
this correlation ([EE]), although with a slightly larger scatter of the data points. We 
have further updated this sample to 40 GRBs with known redshifts and Fig. [El (right 
panel) shows the E pea k-Li SO correlation found with this sample. In Tab.l are reported 
the fit results of the -Speak — L iso correlation. The scatter of the data points around this 
correlation (see Fig lT2|) can be modeled with a Gaussian of er ~ 0.2. Also in this case, 
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as discussed for the Speak — -^iso ([55]) the resulting reduced x 2 is quite large unless 
accounting for a sample variance of the order of ~ 0.13. 

As discussed in |l()8j . the luminosity Lj SO is defined by combining the time- 
integrated spectrum of the burst with its peak flux (also E peak is derived using the 
time integrated spectrum). This assumes that the time-integrated spectrum is also 
representative of the peak spectrum. However, it has been shown that GRBs are 
characterized by a considerable spectral evolution (e.g. 68J). If the peak luminosity is 
derived only by considering the spectrum integrated over a small time interval (~ 1 sec) 
centered around the peak of the burst light curve , we find a larger dispersion of the 
Epeak-Eiso correlation (see |l()8j ). This suggests that, in general, the time averaged 
quantities (i.e. the peak energy of the time integrated spectrum and the "peak- 
averaged" luminosity) are better correlated than the "time-resolved" quantities. 

Schaefer ( |ll()j ) combined the Lag-Luminosity and the Variability-Luminosity 
relations of 9 GRBs to build their Hubble diagram. He showed that GRBs might 
be powerful cosmological tools. However, some notes of caution should be addressed 
when using the correlations presented in Tab.l for cosmographic purposes. In fact, 
one fundamental condition is that there exist at least one cosmology in which these 
correlations give a good fit, i.e. a reduced Xred — 1 an d, unless accounting for the 
sample variance, this is not the case for the correlations presented so far. 

3.5. The isotropic luminosity-peak energy-high signal timescale correlation 

(Lisa --Epeak ~^~0.45 ) 

A new correlation recently discovered by Firmani et al. ([.36 J relates three observables 
of the GRB prompt emission. These are the isotropic luminosity L iso , the rest frame 
peak energy E peilk and the rest frame "high signal" timescale T 045 . The latter is a 
parameter which has been previously used to characterize the GRB variability (e.g. 
[89 ) and represents the time spanned the brightest 45% of the total counts above the 
background. Through the analysis of 19 GRBs, for which L iso -Ep e ak and T0.45 could be 
derived, [SHI found that L iso oc E^ k ■ Tq ® 5 49 with a very small scatter. This correlation 
is presented in FigJHl 

The L iso — -Epeak — T QA5 correlation is based on prompt emission properties only 
and it has some interesting consequences: (i) it represents a new powerful (redshift) 
indicator for GRBs without measured redshifts, which could be computed only from 
the prompt emission data (spectrum and light curve); (ii) it represents a new powerful 
cosmological tool ( |37| - see Sec. 6.2) which is model independent (differently from the 
E 1 — E peak correlation (see Sec. 5) which relies on the standard GRB jet model); (iii) it 
is "Lorentz invariant" for normal fireballs, i.e. when the jet opening angle is 9 jet > 1/r. 
In this case, in fact, the luminosity scales as 5~ 2 when it is transformed from the rest 
to the comoving frame and the doppler factor cancels out with that of the peak energy 
and of T0.45 - see 
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Figure 6. Multi variable linear regression analysis of the ii S o-^peak-?o.45 correlation 
based on 19 GRBs. The open circles represent the 19 bursts used for the fit (solid 
line). The solid filled region represents the la uncertainty on the best fit (see Eq. 1). 
The fit is performed in the centroid defined by the data points (solid circle) where the 
errors on the best fit parameters are uncorrelated. GRB 050525 was not used for the 
multi-variable regression analysis due to its uncertain luminosity. 



4. The third observable: the jet break time (tbreak) 

As shown above the correlations between the spectral peak energy and the isotropic 
energy or luminosity are too much scattered to be used as distance indicators. Only the 
L iso --E peak -To.45, which is based on prompt emission properties, is sufficiently tight to 
standardize GRB energetics (as it will be shown in Sec. 5 and Sec. 6). 

However, all the above correlations have been derived under the hypothesis that 
GRBs are isotropic sources. Indeed, the possibility that GRB fireballs are collimated 
was first proposed for GRB 970508 ([10]) and subsequently invoked for GRB 990123 as 
a possible explanation for their extraordinarily large isotropic energy The main 

prediction of the collimated GRB model is the appearance of an achromatic break in 
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Figure 7. The E pe ak-E^ correlation from obtained under the assumption that the 
circum burst medium has a constant density. The 18 bursts are those with measured 
redshift and firm measurements of their spectral properties. The best powerlaw fit of 
the data points, obtained by accounting for the uncertainties on both coordinates, is 
shown (solid line) together with the linear regression analysis, i.e. without accounting 
for the data uncertainties (dotted line) . The scatter of the data points around the best 
fit correlation is shown in the insert with its Gaussian fit. The 1,2 and 3a scatter of 
the data points around the best correlation is represented by the shaded region. 



the afterglow light curve which, after this break time, declines more steeply than before 
it EH])- Due to the relativistic beaming of the photons emitted by the fireball, the 

observer perceives the photons within a cone with aperture #r oc 1/T, where V is the 
bulk Lorentz factor of the material responsible for the emission. During the afterglow 
phase the fireball is decelerated by the circum burst medium and its bulk Lorentz factor 
decreases, i.e. the beaming angle #r increases with time. A critical time is reached 
when the beaming angle equals the jet opening angle, i.e. 6-p ~ 1/T ~ 0j et , i.e. when 
the entire jet surface is visible . Under such hypothesis the jet opening angle 0j et can 
be estimated through this characteristic time ([HH]), i-e. the so called jet-break time 
^break of the afterglow light curve. Typical tbreak values ranges between 0.5 and 6 days 
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Figure 8. The E voa k-E 1 correlation in the case of a wind profile of the external 
medium density as found with the 18 GRBs reported in The solid line represents 
the best fit powerlaw model obtained accounting for the errors on both coordinates. We 
also show the fit obtained with the simplest linear regression, i.e. without accounting 
for the errors on the coordinates (dotted line), which has a slope of 0.92. The 
circled point represents GRB 030326 which is giving the largest contribution (23%) 
to the best fit reduced x 2 . The shaded regions represent the 1, 2 and 3<r scatter 
around the best fit correlation. The names of the 18 GRBs are indicated. The insert 
reports the distribution (hatched histogram) of the scatter of the data points computed 
perpendicular to the best correlation (solid line in the main plot) and its gaussian fit 
(solid line in the insert) which has a a — 0.08. 

The jet opening angle can be derived from tbreak in two different scenarios (e.g. 
|lllj ): (a) assuming that the circum burst medium is homogeneous (HM) or (b) 
assuming a stratified density profile (WM) produced, for instance, by the wind of the 
GRB progenitor. In both cases the jet is assumed to be uniform. 

5. The peak energy— collimation corrected energy correlation (E pea k-E 7 ) 

By correcting the isotropic energy of GRBs with measured redshifts, under the 
hypothesis of a homogeneous circum-burst medium, Ghirlanda et al. (|34|) discovered 
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Table 2. GRB intrinsic correlations involving collimation corrected quantities. N: 
number of GRBs used to find the correlation. The rank correlation coefficient r s and 
its chance probability P are given. The best fit parameters (with their la uncertainty) 
are also reported together with their x 2 (dof). 



Coll. corrected Correlations 


N 


r B 


P 


<3 


a 


X 2 (dof) 


Log(E pca-k ) = q + aLog(E y , H u) 


18 


0.9 


2.3xl0" s 


-32.36±2.27 


0.69± 0.04 


22.4(16) 


Log(E pcs-k ) = q + aLog(E^ yWM ) 


18 


0.92 


6.8xl0~ 8 


_49.44±3.27 


1.03±0.06 


18(16) 


Log(E peak ) = q + aLog(E iso ) 


18 






-48.05±0.24 


1.93±0.11 


-1.08±0.17 24.2(16) 


+/3£o3(*brcak) 














Log(E pcllk ) = q + aLog{L^. HM ) 


16 


0.9 


4xl0 -4 


-22.8±1.3 


0.51±0.03 


41(14) 


Log(E pea-k ) = q + aLog(L^ yWM ) 


16 


0.9 


7xl0~ 4 


-35.1±2.0 


0.75±0.04 


76(14) 



that the collimation corrected energy E y = E iso (l — cos 9j Ct ) is tightly correlated with the 
rest frame peak energy. With an initial sample of 15 GRBs ([34 ) this correlation results 
steeper ( i.e. -Epeak oc E®- 7 ) than the corresponding correlation involving the isotropic 
energy and with a very small scatter (i.e. a ~ 0.1 to be compared with the ~ 0.2 scatter 
of the correlation with Ei SO ) . The addition of new GRBs (some of which also discovered 
by Swift and, at least, by another IPN satellite in order to have an estimate of the peak 
energy) has confirmed this correlation and its small scatter (e.g. [105] ). 

Recently, Nava et al. ([HE]) reconsidered and updated the original sample of GRBs 
with firm estimate of their redshift, spectral properties and jet break times. With the 
most updated sample of 18 GRBs we can re-compute the E pea ^-E^ correlation in the 
HM case. This is represented in Fig[7|and results -E pea k oc ^°- 69±0 04 ) ( se e Tab. 2) with 
a very small scatter (i.e. a ~ 0.1). Only two bursts (i.e. GRB 980425 and GRB 031203 
- which are not reported in the figures - but see [105] ) are outliers to this correlation, 
as well as to the E peak -E iso correlation. Possible interpretations has been put forward 
such as the possibility that these events are seen out of their jet opening angle ( |112j ) 
and therefore appear de-beamed in both -E pea k and E iso (see also |113| ). However, other 
possible explanations could be considered ( [H4J . ITTti; ) . 

Most intriguingly, [23] also computed the E peak -E^ correlation in the WM case 
where a typical r~ 2 density profile is assumed. The resulting correlation is reported in 
FigGD The -E pea k--E'gamma correlation derived in the WM case has two major properties: 
(i) its scatter, i.e. a = 0.08, is smaller than in the HM case (although, due to the still 
limited number of data points, this does not represent a test for the WM case against 
the HM case) and (ii) it is linear, i.e. -E pea k oc E T Under the hypothesis that our line of 
sight is within the GRB jet aperture angle, the existence of a linear E peg ^-E 1 correlation 
also implies that it is invariant when transforming from the source rest frame to the 
fireball comoving frame. A striking, still not explained, consequence of this property 
is that the total number of photons emitted in different GRBs is similar and should 
correspond to ~ 10 57 i.e. roughly the number of baryons in 1 solar mass. The latter 
property might have important implications for the understanding of the dynamics and 
radiative processes of GRBs. 

One of the major criticism encountered by the above correlations between the peak 
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Figure 9. The 3-D representation of the -E P eak - -E'iso - £break correlation computed with 
the updated sample of The shaded plane represents the best fit correlation and 
the data points are marked differently if the lie above (red points) or below (gray stars) 
this plane. The vertical lines show the distance of the data points from the best fit 
plane. 



energy and the collimation corrected energy is the fact that the collimation correction 
is derived from the measure of tbreak by assuming the standard fireball model. If on the 
one hand the small scatter of the E peak -E 7 correlation might be, in turn, regarded as a 
confirm of this model, on the other hand it could still be a matter of debate when these 
correlations are used for cosmographic purposes. 

6. The peak energy— isotropic energy—jet break time correlation 

( -Speak~-^'iso - ^break ) 

A completely empirical correlation also links the same three observables used in deriving 
the above E peak -E-y correlation. Liang & Zhang ([2H1), in fact, found that the peak 
energy -E pea k, the jet break time tbreak and the isotropic energy i? iso are correlated. 
In FigEl we report the 3-D correlation obtained with the most updated sample of 18 
GRBs for which all the three observables are firmly measured. The scatter of the data 
points around this correlation allows its use to constrain the cosmological parameters 
([3U]). As discussed in [33] the model dependent Ep^-E^ correlations (i.e. derived 
under the assumption of a standard uniform jet model and either for a uniform or a 
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wind circumburst medium) are consistent with this completely empirical 3-D correlation. 
This result, therefore, reinforces the validity of the scenario within which they have been 
derived, i.e. a relativistically fireball with a uniform jet geometry which is decelerated 
by the external medium, with either a constant density or with an r~ 2 profile. 

Similarly to what has been done with the isotropic quantities, we can explore if 
the collimation corrected E pea k-E^ correlation still holds when the luminosity, instead 
of the energy, is considered. In the second part of Tab. 2 we report also the correlations 
between -Speak and L 7 in the HM case (which represents an update of the correlation 
discussed in |lU8j ) and also the new correlation in the WM case. We note that both the 
HM and the WM correlations defined with L 7 are more scattered than the corresponding 
correlations defined with the collimation corrected energy E^. 

6.1. Probing the E pea ^ — E [so correlation 

One of the major drawback of the -E pea k — -Eiso correlation is that it has been found 
with a small number of GRBs. This led to think that some selection effects might play 
a relevant role in deriving this correlation. In support of this correlation we should 
mention that the extension ( [33] , jSS] ) of the original sample of GRBs from which it was 
derived (|2>J) leads to find again a statistically significant correlation between -Speak and 
E- 

-'-'ISO- 

One test that can be performed even without knowing the redshift of GRBs is to 
check if bursts with measured fluence S and peak energy -E°^ k are consistent with the 
-Epeak — E iso correlation, by considering all possible redshifts. The application of this test 
to a sample of BATSE bursts ( |103| . |l()4j ) led to claim that the E pea k-Ei SO correlation 
is not satisfied by most of the BATSE GRBs without a redshift measurement. Instead, 
a different conclusion has been reached by [107J. 

We performed (108J) a different test by checking if a large sample of GRBs 
with a pseudo-redshift measurement could indicate the existence of a relation in the 
Epeak — E [so plane, although with a possible larger scatter than that found for the E pe3 k- 
_E iso correlation with few GRBs with spectroscopic measured redshifts. 

With a sample of 442 long duration GRBs whose spectral properties has been 
studied and whose redshifts has been derived from the spectral lag-luminosity 

correlation ( |106j ). we populated the -E pea k — E J;iso plane (black crosses in Fig. [TUJ). We 
found that this large sample of burst produce an E pea k-E iso correlation (solid line in 
FiglTUI) whose slope (normalization) is slightly flatter (smaller) than that found with 
the sample of 43 GRBs (red filled circles in FigHnj). The gaussian scatter of the pseudo- 
redshift bursts around their correlation has a standard deviation a = 0.2, i.e. consistent 
with that of the 43 GRBs around their best fit correlation (long-dashed red line in 
FigHnj). This suggests that indeed a correlation between -Epeak and -E iso exists. 

However, it might still be argued that the correlation found with the 43 and with 
the 442 GRBs have inconsistent normalization, although similar scatter and slopes. This 
is shown by the 43 GRBs (red filled circles in FiglTUI) with spectroscopically measured 
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Figure 10. Rest frame peak energy -Epeak versus isotropic equivalent energy E- lso . 
Red filled circles are the 43 GRBs with spectroscopically measured redshifts and 
published spectral properties (adapted from |SS|)- The long-dashed line is their best 
fit (weighting for the errors on both variables) which is reported in Tab. 1. The 
black crosses are the 442 GRBs with pseudo redshifts derived from the lag-luminosity 
relation. The solid line is the best fit to these data points which results -E pe ak oc 7 
with a reduced Xi = 4-0 (440 dof). The shaded region represents the 3a scatter of 
the black points around their best fit line (solid line). The triple-dot-dashed line 
represents the 2? pca k — correlation. 



redshifts which lie on the right tail of the scatter distribution of the 442 GRBs (black 
crosses) in the -E pea k — E [so plane. This fact might be interpreted as due to the selection 
effect of the jet opening angle. In fact, from the scatter of the 442 GRBs we can derive 
their jet opening angles by assuming that the Ep^-E^ correlation exists and that its 
scatter is (as shown in [HI]) much smaller than that of the E pe3 k — -^iso correlation. The 
8 j et distribution derived in this way is shown in Fig^^ (solid histogram) and it is well 
represented by a log-normal distribution (solid line) with a typical 9j et ~ 6°. The 9- ]et 
distribution of the 19 GRBs with measured tj et is also shown (dot-dashed histogram 
in FigHTJ) and it appears shifted to the small-angle tail of the 9j et distribution of the 
442 angles. This suggests that the 43 GRBs, which are used to define the -E pea k — E iso 
correlation, have jet opening angles which are systematically smaller than the average 
value of the distribution and this is what makes them more luminous and, therefore, 
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Figure 11. Jet opening angle distributions. The solid histogram represents the #j e t 
derived from the large sample of 442 GRBs with pseudo redshifts requiring that they 
satisfy the -Epeak — E~ correlation ( |l(J5) .p?5 ] - with the approximation of no scatter). 
The solid line is the best fit log-normal distribution. The dot-dashed histogram 
represents the angle distribution of the 19 GRBs with spectroscopic redshifts and 
well constrained ij e t- 

more easily detected. 

Interestingly the fact that the scatter of the 442 GRBs in the £^p ea k — -^iso plane 
around their best fit correlation is gaussian indicates that, if the scatter is due solely 
to their angle distribution, this is characterized by a preferential angle. If, instead, the 
angle distribution were flat, we should expect a uniform distribution of the data points 
around the -Ep ea k — -^iso correlation. 

7. Constraining the cosmological parameters with GRBs 

The scatter of the correlations described in the previous sections are compared in FigIT2l 
The collimation corrected correlations E pea ] s -E li HM and E pe3 k-E J} wM or the empirical 
correlation -Epeak^-Eiso^break, due to their very low scatter can be used to constrain the 
cosmological parameters. Moreover, all these correlations have an acceptable reduced 
X 2 in the concordance cosmological model. 

As a first test of the possible use of the above correlations for cosmography we 
show in Fig. EH the luminosity distance derived with the Ep^—E^ correlation in the 
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Figure 12. Left: Distribution of the scatter around the correlations involving 
isotropic quantities (reported in Tab.l). The orange distribution (and solid orange 
line) represents the scatter distribution (and its gaussian fit) around the -E p0 ak- 
i?i so correlation. The blue histogram and violet histograms are for the -Speak - 
L; so and V — Li ao correlation, respectively. Right: Distribution of the scatter 
around the correlations involving the collimation corrected quantities and of the 3- 
D empirical correlation (reported in Tab. 2). The green distribution (and dashed green 
line) represents the scatter distribution (and its gaussian fit) around the E pea k-E~f 
correlation (in the wind circumburst scenario). The orange and pink histograms are 
for the scatter of the Ep^-Ey (in the homogeneous circumburst medium scenario) 
and E lso — E pca _k — ibreak correlation, respectively. 

HM and WM scenario for the 18 GRBs of the sample discussed in [32]. By comparing 
the plots with Fig. El we note that by using the E- pea \ r E 1 correlation the dispersion of 
the data points is reduced with respect to the assumption that all GRBs have a unique 
luminosity (see Sec. 2). Also note that no substantial difference appears if using the 
Epeak-Ery correlation derived either in the HM or in the WM case (top and lower panel 
of Fig. IT31 respectively). 

Different methods have been proposed to constrain the cosmological parameters 
with GRBs (001) |H6| . |41j). A common feature of these methods is that they construct 
a x 2 statistic, based on one of the correlations described before, as a function of two 
cosmological parameters. The minimum \ 2 an d the statistical confidence intervals (i.e. 
68%, 90% and 99%) represent the constraints on the couple of cosmological parameters. 

One important point (see also Sec 7.3) is that the correlations that we have 
presented so far have been derived assuming a particular cosmological model, i.e. 
Om=0.3, ^a=0.7 and h=0.7. Indeed, in the E pe3 k-E J} HM correlation, for instance, 
-Speak = -Speak(l + 2 i s independent from the cosmological model whereas depends on 
the cosmological parameters through the luminosity distance, i.e. = f (Att D^E) / (1 + 
z), where D L is a function of (0 M ,OA,/t) (see Sec. 1.1). The need to assume a set of 
cosmological parameters to compute the above correlations is mainly due to the lack 



Cosmology with GRBs 



22 



29. 5F 




z 



Figure 13. Hubble diagram of the 18 GRBs of the sample of derived with 
the Epeak-Ej correlations in the HM and WM case. The curves represent different 
cosmological models: (f2^f > ^a)=(0.3,0.7) [solid-blue line], (fi^f, fiA.)=(0. 0,1.0) [long- 
dashed red line], (Om, ^a) = (1.-0,0.0) [short-dashed green line]. 

of a set of low redshift GRBs which, being cosmology independent, would calibrate 
these correlations. Therefore, if we want to fit any set of cosmological parameters with 
the above correlation we cannot use a unique correlation i.e. obtained in a particular 
cosmology. This method, which has indeed been adopted by |116j . would clearly suffer 
of such a circular argument. 

The "circularity problem" ( |117j ) can be avoided in two ways: (i) through the 
calibration of these correlations by several low redshift GRBs (in fact, at z < 0.1 the 
luminosity distance has a negligible dependence from any choice of the cosmological 
parameters Qm^a) or (h) through a solid physical interpretation of these correlation 
which would fix their slope independently from cosmology. A third possibility, presented 
in Sec. 7, consists in calibrating the correlations with a set of GRBs with a low dispersion 
in redshift. 

While waiting for a sample of calibrators or for a definitive interpretation, we can 
search for suitable statistical methods to use these correlations as distance indicators 
but avoiding the "circularity problem". The methods that have been used to fit the 
cosmological parameters through GRBs can be summarized as follows: 

• (I) The scatter method ([ID]). By fitting the correlation for every choice of the 
cosmological parameters that we want to constrain (e.g. Qm, Ra.) a x 2 surface (as 
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a function of these parameters) is built. The best cosmological model corresponds 
to the minimum scatter around the correlation (which has to be re-computed in 
every cosmology in order to avoid the circularity problem) and is identified by the 
minimum of the x 2 (^m,^a) surface. 

• (II) The luminosity distance method ( [4*0] . [TTT)] . [3H] ) . The main steps are: (1) choose 
a cosmology and fit the E pea \ s -E 1 correlation ; (2) from the best fit correlation 
the term E 1 is estimated; (3) from the definition of E 1 = / • Ei SO (where also 
/ = 1 — cosset is a function of E iso through 6- ]et ) derive E iso from which (4) the 
luminosity distance Dl, c is computed; (5) build a x 2 by comparing Dl^ c with that 
derived from the cosmological model D^. By repeating these steps for every choice 
of the cosmological parameters a \ 2 iS^M , ^a) surface is derived. Also in this case the 
minimum x 2 represents the best cosmology, i.e. the model which best matches the 
luminosity distance derived from the correlation and that predicted by the model 
itself. 

• (III) The Bayesian method (|H]). The two methods described above are based 
on the concept that a correlation exists between two variables, e.g. -E pea k and E 1 . 
However these correlations are also very likely related to the physics of GRBs and, 
therefore, they should be unique. This property is not exploited by methods (I) and 
(II). Firmani et al. (|41j) proposed a more sophisticated method which makes use 
of both these conditions, i.e. the existence and uniqueness of the correlation. The 
method is based on the Bayes theorem (e.g. |118| ). The basic steps are: (1) choose 
a cosmology Cl and fit the correlation; (2) test this correlation (i.e. keep it fixed) in 
all the possible cosmologies Q and derive a conditioned probability surface P(Q\Q) 
which represents a "weight" for the starting cosmology Cl (i.e. it describes how 
all the cosmologies that have been tested "respond" to the starting cosmology); 
(4) by repeating these steps for different starting cosmologies Q and combining, 
at each step, the new conditioned probability with those found in the previous 
iterations, the Bayesian methods naturally converges to a final probability surface 
whose maximum represents the best cosmological model. 

The last method, which exploits all the information contained in the adopted 
correlation, has also a practical advantage: it solves the problem of the loitering line 
which separates Big-Bang and no-Big-Bang universes. As shown in Fig. El near this line 
the luminosity distance stripes are highly wounded-up and clustered. In correspondence 
to this line in fact Dl falls rapidly to zero. If the adopted sample of standard candles 
extends to very high redshifts, this region can attract the minimum of the x 2 surface, as 
shown in Fig.l of [HI]. For this reason the Bayesian method is more accurate because 
it assigns a low probability to the cosmologies near this line (see |41j). 

8. Constraints on Qm and Q\ 

With the Bayesian method we can constrain the cosmological parameters Qm and Q\ 
with the sample of 18 GRBs reported in ^ updated with GRB 051022 (see Eg). 
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Figure 14. Regions of ±1% variation around lines of constant d-^ in the (Om^a) 
plane at redshift 0.5,1,1.5 and 2. assuming that each line passes through the fiducial 
pomt(n M =0.33,fi A =0.77). 



In Fig. HSlwe show the contours corresponding to the 68%, 90% and 99% confidence 
levels obtained with all the three correlations, i.e. E peak -E 7 ^ HM (solid line), -Ep Ca k- 
E Jt wM (dashed line) and the empirical 3D correlation -Epeak^-Eiso^break (dotted line). 
We note that the contours obtained with these correlations (of which the first two are 
model dependent and the third is completely empirical) are consistent one with another. 
We also note that the centroid of the confidence contours are only slightly different. 
However, given the large allowance area of the same contours the exact position of the 
centroid is not relevant, i.e. the contours obtained with GRBs as standard candles are 
consistent with the concordance model. 

8.1. Statistical errors 

One of the critical issue in the use of the three correlations for cosmology is related 
to statistical errors. There are, in fact, two major sources of errors which shape the 
contours obtained with GRBs: the statistical errors associated with the observables 
used to calculate E peilk and E 1 and the errors on the best fit parameters (slope and 
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Figure 15. Constraints on the Q,m and £1a parameters obtained with the Bayesian 
method through the three correlations: Ep^-E^HM (solid red line), E pea ^-E^ wM 
(dashed green line) and the empirical correlation Ep Ca ^-E- lso -t' hrca ^ (dotted blue line). 
The lines correspond to the 68%, 90% and 99% confidence contours. 

normalization) of the correlation resulting from the fitting procedure. These errors 
enters in the construction of the x 2 surface. 

In particular, the issue of accounting for all the sources of errors makes the contours 
particularly large as in the case of the empirical 3D correlation which involves three 
observables (with their statistical uncertainties) and three fit parameters (with their 
errors, resulting from the best fit of the data points). As a simple test we compare 
in Fig. the constraints on Qm and when accounting only for the errors on the 
observables (solid-red contours) or when accounting also for the errors on the best 
fit parameters (dashed blue contours). In Fig. (left panel) we show the constraints 
obtained with the empirical E pea k-E iso -threak correlation by ignoring (solid contours) the 
errors associated to the the best fit parameters which in this case are the normalization 
and 2 slopes (because the correlation is expressed, for instance, as Ei SO as a function of 
.Speak and ibreak)- The same figure shows, instead, that if we account for all the errors, 
we find (dashed blue contours) considerably larger constraints. This difference is also 
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represented in the right panel of Fig. EH where the error on the distance modulus \x 
obtained accounting for all the errors is represented versus the same quantity obtained 
including only the errors on the observables. 
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Figure 16. Left: Constraints on the cosmological parameters Qm,Qa obtained with 
the empirical 3D correlation and the updated sample of 19 GRBs. The solid (red) 
contours are obtained by considering, as in |39| . only the errors on the variables, i.e. 
-Epeak, -Eiso and ^brcak- The long-dashed blue contours include also the errors on the 
best fit parameters. The dotted contours are obtained with the original sample of 14 
GRBs used by [21], accounting only for the errors on the variables. Right: Errors 
on the distance modulus /i computed considering the errors on the best fit correlation 
parameters (k,A) versus the same error computed ignoring these uncertainties and 
considering only the errors on -E 7j i SO , E p and <brcak- In this case we assumed that the 
exponent of tbreak is linear and we ignored the error on this parameter. 

The dotted (red) contours in Fig. EH (left panel) are obtained with the sample of 
14 GRBs reported in the original work of [39J without including the errors on the best 
fit parameters. Here we adopted the Bayesian method to derive all the cosmological 
contours. Moreover, another possible reason for the difference between our contours 
(red dotted in Fig. HT)|) and the slightly smaller contours reported in [HH] might be the 
fact that our method for fitting the 3D empirical correlation takes into account the 
errors on the three variables and adopts a the \ 2 minimization method to find the best 
fit parameter values, whereas the multiregression analysis ([20]) does not account for all 
the errors on the observables. 

8.2. Combining GRBs and SNIa 

One important consequence of using GRBs as standard candles through the E peg ^- 
E 1 correlations or the empirical 3-D £ , iso -£ , pea ] (; -tbreak correlation is that they can be 
combined with SNIa. With respect to SNIa, GRBs already extend in a redshift range 
which is beyond the present (and future) SNIa redshift limit (i.e. z ~ 1.7 predicted 
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Figure 17. Constraints on the cosmological parameters Qm, obtained with the 
updated sample of 19 GRBs, in the homogeneous density case (HM). The solid (red) 
contours, obtained with the 19 GRBs alone, represent the 68.3%, 90% and 99% 
confidence regions. The center of these contours (red cross) corresponds to a minimum 
X 2 = 15.25/17 dof and has £Im = 0.23 and £Ia = 0.81. The contours obtained with 
the 156 SN la of the "Gold" sample of are shown by the dashed (blue) lines. The 
joint GRB+SN constraints are represented by the shaded contours. We also show the 
90% confidence contours obtained with the WMAP data (from 6 ). 

for SNAP - e.g. [119J). Furthermore, GRBs are detected in the 7 ray band which is 
unaffected by dust extinction limitations. However, SNIa are detected in large number 
also at very low redshifts and their "stretching-luminosity" correlation can be well 
calibrated. Therefore, GRBs should be considered as complementary to SNIa at very 
high redshisfts. 

In Fig. El we show the contours obtained with GRBs (in the HM case) and we 
combine them with the constraints obtained for a large sample of 156 SNIa (the "gold" 
sample of JI]). The combined fit is clearly dominated by the large number of SNIa 
with respect to GRBs (a factor 10 more numerous), however, the power of GRBs is to 
make the joint contours more consistent with the concordance model, also due to the 
different orientation of the contours obtained by these two probes (see Sec. 1). 
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Figure 18. Contours at 68.3%, 95.5%, and 99.7% CL's obtained from projecting 
to the rijvi , f^A plane the iiscr-Epeak-^o^s relation that is obtained from the fit of 
the GRB data at each value of the (Om^a) pair. This plot shows that the L; so - 
-E"peak-2b.45 relation is sensitive to cosmology, so that it may be used to discriminate 
cosmological parameters if an optimal method to circumvent the circularity problem is 
used. The diagonal line corresponds to the flat geometry cosmology, the upper curve 
is the loitering limit between Big Bang and no Big Bang models, and the lower curve 
indicates the division between accelerating and non-accelerating universes. 

8.3. Cosmological constraints with the Li SO -E pea k-ToA5 correlation 

Finally, we present in Fig. [T^l the constraints ( 37 ) on the parameters f2 M and 
obtained with the "prompt-emission" correlation L iso -.Ep ea k-To.45 (see Sec. 3.5). 

First we can apply the scatter method (I) to test if this correlation is, indeed, 
sensitive to the cosmological parameters. We indeed show in Fig. EH smallest scatter is 
obtained for (f2M,^A) = (0.3,0.7), i.e. very close to the concordance model. 

By adopting the same Bayesian method described above, we note that the 
constraints obtained with this new correlation with a sample of only 19 GRBs are 
tighter than those obtained with the E peak -E^ correlation fFig. fTTjl . 

The best cosmological model, indicated by the star symbol in Fig. EH corresponds 
to flyi = 0.3lio;gg,fiA = 0.80^o;3Q (lcr uncertainty), i.e. vary close to the flat geometry 
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Figure 19. Constraints on the (Om^a) plane for a A cosmology from the GRB 
Hubble diagram using our Bayesian method to circumvent the circularity problem and 
from the gold set of SNIa of (thin solid lines). The GRB constraints are obtained 
with the £iso-£pcak-7o.45 correlation. 

case. Indeed, if a flat geometry is assumed, we find f^M = 0.29to;og (P^P- 
9. The future of GRB-cosmology 

The luminosity distance test requires large samples of sources possibly distributed over 
a wide redshift range. This can be understood from Fig. E] which shows the topology 
of the luminosity distance as a function of the cosmological parameters that we want 
to constrain. In this case we consider Qm arid Q\ as the two free parameters. The 
stripes in Fig. El show the degeneracy of the Qm and Qa parameters and how this 
degeneracy varies with the redshift z. Each stripe represents, for a fixed redshift z, all 
the possible combinations (Q M ,Q^) which give a luminosity distance equal within 10%. 
By increasing the redshift the stripes wound up due to the topology of the luminosity 
distance as a function of these two parameters, i.e. (Om 5 Ra.) ( see |H7| ). It is clear from 
Fig. E] that if we use a sample of sources distributed over a wide range of redshifts 
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we "intersects" differently oriented stripes and we end up with accurate constraints on 
the two cosmological parameters. On the other hand a sample of sources distributed 
in a very small redshift range would not break the degeneracy of these two parameters. 
The requirement of a large number of sources instead is a key ingredient to reduce the 
possible effect of systematic errors. 

The present open issues related to the use of GRBs in cosmologies are: (1) the low 
number of events: the sample of GRBs which can be used to constrain the cosmological 
parameters through the above correlations are in fact 19 (to date); (2) the lack of a 
calibration sample, i.e. the absence of GRBs at very low redshifts in the sample used 
for cosmology, does not allow to calibrate the correlation and requires to adopt a method 
to fit the cosmological parameters which avoids the circularity problem; (3) the possible 
lensing effects which might introduce a dispersion of the luminosities and energetics 
of GRBs at very high redshifts. The latter point, however, if present, could be easily 
recognized: if strong lensing by compact objects (if by clusters they should be observed) 
is present, we should observe a repetition of a similar structure (e.g. peak shape) within 
the same light curve. Moreover, the lack of a theoretical interpretation of the physical 
nature of all these correlations represents a still open issue. 

The increase of the number of bursts which can be used to measure the cosmological 
parameters, and the possible calibration of the correlations would greatly improve the 
constraints that are presently obtained with few events and with a non-calibrated 
correlation. 

In order to use GRBs as a cosmological tools, through the above correlations, 
three fundamental parameters, i.e. -E pea k, E 1 and tj et , should be accurately measured. 
This requirement also applies to the case of the empirical correlation of jHSj- On the 
other hand the L iso -E peak -T QA5 does not require the knowledge of the afterglow emission 
because it completely relies on the prompt emission observables. 

However, only a limited number of GRBs, i.e. 19 out of ~ 70 with measured z, 
can be used as standard candles. Presently the most critical observable is the peak 
energy -E pea k- In fact, the limited energy range (15-150 keV) of BAT on board Swift 
allows to constrain with only moderate accuracy the -E pea k of particularly bright-soft 
bursts. However, given the perspective of the cosmological investigation through GRBs, 
it is worth exploring the power of using GRBs as cosmological probes. The other still 
open issue related to the use of GRBs as standard candles is the so called "circularity 
problem" (see Sec. 5). 

To the aim of calibrating this spectral-energy correlation, GRBs at low redshift 
are required. In fact for z < 0.1 the difference in the luminosity distance computed 
for different choices of the cosmological models [for f2 M ,^A G (0,1)] is less than 8%. 
However, if (long) GRBs are produced by the death of massive stars, they are expected 
to roughly follow the cosmic star formation history (SFR) and we should expect that 
the rate of low redshift events is quite small at z < 0.1. Indeed, a considerable number 
of GRBs with z > 1 should be collected in the next years by presently flying instruments 
(Swift and Hete-IT). At such large redshifts, the cosmological models starts to play an 
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Figure 20. Distributions of redshift z, jet break time tb, observed peak energy -Ep^k 
and isotropic equivalent energy E iso for the 150 simulated bursts (hatched histograms). 
Also show (solid filled histogram) are the distributions of the sample of 19 GRBs used 
to derive the £p C ak — E iso correlation. 



important role. However, if it will be possible to have a sufficient number of GRBs with 
a similar redshift, it might still be possible to calibrate the slope of these correlations 
also with high redshift GRBs (see Sec. 7.3). 

9.1. The simulation 

In order to fully appreciate the potential use of GRBs for the cosmological investigation, 
we simulate a sample of bursts comparable in number to the "Gold" sample of 156 SN 
la. Similar simulations have been presented before (e.g. [121,1321 )• However, different 
assumptions can be made on the properties of the simulated sample and the results 
are clearly dependent on these assumptions. In particular, simulations based on the 
observed parameters (|42j) strongly depend on the almost unknown selection effects 
which shape the observed distributions. 

We adopt here a method which makes use of the intrinsic properties of GRBs as 
described by the -E pea k--E 7 and the £ l P eak--E'iso correlations. We use the most updated 
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version of these correlations as found with the sample of 19 GRBs described in the 
previous sections. 

The assumptions of our simulation are: 

• we assume that GRBs have an "isotropic energy" function which is described by 
a powerlaw N(E iso ) oc E? so for Ei so ,min < E iso < -Eiso.max (e.g. |HH] ) - Further we 
assume that GRBs follow the cosmic star formation rate (SFR - as modeled in 

D2DI ); 

• we use the E peak -E iso correlation to derive the peak energy -E pea k and we model the 
scatter of the simulated GRBs around the E pea k-Ei SO correlation with a gaussian 
distribution with cr = 0.2 (which corresponds to the present scatter of the 19 GRBs 
around their best fit correlation); 

• we use the Ep eak -Ey correlation as found with the 19 GRBs in the WM case to 
calculate Ej and we model the scatter of the simulated GRBs around the Ep^-E^ 
correlation with a gaussian distribution with o = 0.08 (as found in Sec. 4); 

• we derive the jet opening angle 0j and the corresponding jet break time tbreak! 

• we assume that the simulated GRB spectra are described by a Band model spectrum 
( [12 lj ) with typical low and high energy spectral photon indices a = —1.0 and 
j3 = —2.5 and require that the simulated GRB fluence in the 2-400 keV energy 
band is above a typical instrumental detection threshold S ~ 10~ 7 erg/cm 2 . This 
corresponds roughly to the present threshold of Hete-II in the same energy band. 

Following the procedure described above we built a sample of 150 GRBs with the 
relevant parameters: z, E- mo , E peak , tbreak- The errors associated to these parameters 
are assumed to be cosmology invariant and they are set to 10%, 20% and 20% for 
-Else -Epeak, ^break respectively. Moreover, we model the GRB intrinsic isotropic energy 
function with a powerlaw with 5 = —1.3 between two limiting energies (10 49 -10 55 erg). 
This particular choice of parameters is due to the requirement that the distributions of 
the relevant quantities (shown in Fig. l2*Uj) of the simulated sample are consistent with 
the same distributions for the present sample of 19 GRBs (red histograms in Fig. EDI) - 
The sample is simulated in the standard cosmology (Qm — 0.3 and fl\ = h = 0.7). 
We compare the distributions of z, U, E p eak and E iso for the 150 simulated bursts with 
the same distributions of the 19 GRBs in Fig. |20J We note that by choosing a steeper 
energy function we obtain a much larger number of XRF and XRR with respect to 
normal GRBs. 

9.2. Cosmological constraints with the simulated sample 

The results obtained with the sample of 150 simulated GRBs are presented in Fig. IIH] in 
this case the constraints are comparable with those obtained with SNIa. By comparing 
the la contours of GRB alone from Figl^TJ to the same contours (solid line) of Fig. El 
(obtained with the 19 GRBs), we note an improvement (of roughly a factor 10) with 
the sample of 150 bursts. 
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Fig. |^ also shows the different orientation of the GRB contours with respect to 
SN la due to the "topology" of the luminosity distance as a function of the Qm-^a 
parameters. Most GRBs of our simulated sample z ~ 2, and this 

explains the tilt of the contours in the Q plane. Clearly the contours obtained through 
the simulated sample depends on the assumptions: in particular we have no knowledge 
of the burst intrinsic energy function N(E iso ). However, if we accept the hypothesis to 
model it with a simple powerlaw, we can change the slope and also include the possible 
effect of the redshift evolution. We tested the dependence from these assumptions of 
the constraints reported in Fig. |^and found that, by assuming different S values and 
a (1 + z) evolutionary factor, what changes is the redshift distribution of the simulated 
sample and therefore the orientation of the GRB contours in Fig. |^ The same happens 
if we adopt, for the same choice of parameters reported above, a different SFR. 

Further, we can also use the CMB priors. First we assume the cosmological constant 
model with the 2 CMB priors, i.e (i) Vt tot = 1 and (ii) Qm — 0.14/h 2 . In this case the only 
free parameter is h (or equivalently f2ju) We obtain the best fit values of Qm = 0.27±0.02 
and Q A = 0.73 ± 0.02. 

One of the major promises of the cosmological use of GRBs is related to the 
possibility to study the nature of Dark Energy with such a class of "standard candles" 
extending out to very large redshifts. With the present sample of 19 GRBs we can 
explore the equation of state (EOS) of DE, which can be parametrized in different 
ways. Given the already considerably large dispersion of GRB redshifts (i.e. between 
0.168 to 3.2 for the 19 GRBs of our sample) we adopt the parametrization proposed by 
[T22] for the EOS of DE, i.e. P = w(z)p, where: 



With this assumption the luminosity distance, as derived from the Friedman equations, 
is 



which depends on the (wo,w a ) parameters. Note that Eq. El is derived with the prior of 
a flat Universe. 

First, we can assume the CMB prior of a flat universe together with the assumption 
of a non-evolving equation of state of the Dark Energy (i.e. w a = 0). We show in Fig. 
5 the contours obtained with the sample of 150 simulated GRBs and compare with the 
same constraints derived with the 156 SN la of the "Gold" sample. The constraints on 
w a ,w are reported in Fig. 6, assuming Q M = 0.3. 

9. 3. The calibration of the correlations 

The cosmological use of the E peak = K ■ E-* correlation suffers from the so called 
"circularity problem" : this means that both the slope g and the normalization K 



w(z) = w + 



l + z 



(1) 
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Figure 21. Wind density profile case. Constraints on the cosmological parameters 
f^Mj^A obtained with the sample of 150 GRBs simulated by assuming the Ep^-E^ 
correlations derived in the wind density profile. The solid (red) contours, obtained 
with the 19 GRBs alone, represent the 68.3%, 90% and 99% confidence regions on 
the pair of cosmological parameters. The contours obtained with the 156 SN la of 
the "Gold" sample of [TT] are shown by the dashed (blue) lines. The joint GRB+SN 
constraints are represented by the shaded contours. 

of the correlation are cosmology dependent. 

In principle this issue could be solved (a) by a large sample of calibrators, i.e. low 
redshift GRBs for which the luminosity distance Dl is practically independent from 
the cosmological parameters, or (b) by a convincing theoretical interpretation of the 
physical nature of this correlation. In both cases the slope of the correlation would be 
fixed. 

Case (a) could be realized with 5-6 GRBs at z < 0.1. However, if (long) GRBs 
are produced by the core-collapse of massive stars, their rate is mainly regulated by 
the cosmic SFR and, therefore, the probability of detecting events at z < 0.1 is small 
(~ 2 x 10~ 5 ). This number should be convolved with the GRB luminosity function: with 
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Figure 22. Left: Constraints on the cosmological parameters Wq, w a obtained with 
the 150 simulated GRBs (red contours) compared with the same contours obtained 
with the 156 SNIa of the "Gold" sample. A flat universe is assumed (O to t = !)• 
Right: Constraints on the cosmological parameters Wo, obtained with the 150 
simulated GRBs (red contours) compared with the same contours obtained with the 
156 SNIa of the "Gold" sample. A flat universe is assumed (Otot = !)■ 



the assumptions of our simulation we estimate that ~ 1.3% of the 150 GRBs should be 
at low redshifts (i.e. z < 0.4). Instead, we should expect to have more chances to detect 
a considerable number (up to ~ 31%) of intermediate redshift GRBs {z ~ 1 — 2) where 
the cosmic SFR peaks. 

For this reason we explore the possibility to calibrate the correlation using a 
sufficient number of GRBs within a small redshift bin centered around any redshift. 
In fact, if we could have a sample of GRBs all at the same redshift the slope of the 
Epeak-E-f correlation would be cosmology independent. Our objective is, therefore, to 
estimate the minimum number of GRBs (N) within a redshift bin (dz) centered around 
a certain redshift (z c ) which are required to calibrate the correlation. 

In practice the method consists in fitting the correlation for every choice of Q using 
a set of N GRBs distributed in the interval dz (centered around z c ). We consider the 
correlation to be calibrated (i.e. its slope to be cosmology independent) if the change 
of the slope g is less than 1%. 

The free parameters of this test are the number of GRBs N, the "redshift slice" dz 
and the central value of the redshift distribution z c . By Monte Carlo technique we use 
the sample simulated in Sec. 7.1 under the WM assumption to minimize the variation 
of Ag(Q; N, dz, z c ) over the fi M , Vt^ G (0, 1.5) plane as a function of the free parameters 
(N, dz,z c ). 

We tested different values of z c and different redshift dispersions dz G (0.05,0.5). 
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Figure 23. Calibration of the E pea k-E9 correlation. For different samples of GRBs (6, 
12, 18 - corresponding to the solid, dotted and dashed lines respectively) we show the 
maximum variation Ag of the slope of the correlation for any cosmology Q € (0, 1.5) 
as a function of the redshift dispersion of the GRBs dz. The dot dashed line represents 
the limit of variation of 1% of the slope of the correlation. Data points have been 
shifted along the abscissa for graphical purposes. 



We required a minimum number of 6 GRBs to fit the correlation in order to have at least 
4 degrees of freedom. We report our results in Fig. 1221 We show the variation of Ag as 
a function of dz for different samples (6, 12, 18 GRBs - solid, dotted and dashed curves 
in Fig. I2*3*j) . The error bars show the width of the distribution of the simulation results 
and not the uncertainty on the average value. At any redshift the fewer the number 
iV of GRBs the larger the change of Ag (for the same dz) because the correlation is 
less constrained. The dependence from z c is instead different: for larger z c we require a 
smaller bin dz to keep Ag small. 

From the curves reported in Fig. |2H1 we can conclude that already 12 GRBs with 
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z G (0.9,1.1) might be used to calibrate the slope of the E peak -E 7 correlation. At 
redshift z c = 2 instead we require a smaller redshift bin i.e. z G (1.95, 1.05). We find 
that N = 12 GRBs with z G (0.45, 0.75) can be used to achieve the same 1% precision 
in the calibration. However, one key ingredient is that the GRBs used to calibrate the 
correlation do not have the same peak energy otherwise they would collapse in one point 
in the E pea ^-E 1 plane. Within the present sample of 19 GRBs there are only 4 GRBs 
within the redshift interval 0.4-0.8 (i.e. 050525, 041006, 020405 and 051022) and 2 of 
these (050525 and 041006) have a very similar E pea ^. 

We conclude stressing that the three correlations that can be used to constrain 
the cosmological parameters, i.e. the Ep^-E^ (either in the HM and WM case - 
[HH Ell), the empirical -Ei S o~-Epeak~£break correlation (J3S|) and the L iso -i?p eak - 
Tqa5 correlation (jHEllSZl), all involve the peak energy E'peak of the GRB prompt emission 
spectrum. This observable can be properly measured with a detector operating over a 
wide energy range extending up to few MeV as it could be conceived by future missions 
dedicated to collect and study the prompt and afterglow properties of GRB to be used 
as standard candles for cosmology (|HI|). 
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